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^ ■ Abstract 

[ Color diffusion is shown to be an important dissipative property of quark-gluon 

plasmas that rapidly damps collective color modes. We derive the characteristic color 
ff^ , relaxation time scale, tc ~ {SusT log{mE /rnM))~^ , showing its sensitivity to the ratio 

I of the static color electric and magnetic screening masses. This leads to a surprisingly 

small color conductivity, ctc ~ '^T / \og{mE / tum) ^ which in fact vanishes in the semi- 
classical (1-loop) limit. 

Ij 

The quark-gluon plasma (QGP) phase of QCD matter has been predicted to exhibit 
^ \ collective behavior similar to QED plasmas with Debye screened (color) electric interactions, 

^ \ collective plasmon oscillations, and collective color mode instabilities under a variety of 

non-equilibrium configurations (see, e.g., reviews in 0, H). Such collective behavior of 
quark-gluon plasmas (QGP) may be observed experimentally via ultra-relativistic heavy- 
ion reactions if the time scales of dissipative processes that act to damp collectivity are 
sufficiently long. Thusfar, most studies of transport properties of QGP have focused on 
momentum relaxation processes to estimate viscosity and thermal conductivity coefficients, 
dE/dx, etc. 0, ||, |]. Relaxation phenomena associated with the color degrees of freedom 
have received, on the other hand, less attention. In this letter, we derive a new transport 
coefficient, called color diffusion, that has no abelian analog. The related color relaxation 
time, tc, measures how rapidly the color of a parton precesses in an ambient fluctuating 
background field. 
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The relevance of can be seen most clearly in the context of color conductivity. In the 
relaxation time approximation, classical non-abelian transport theory leads to the following 
formula for the static color conductivity 0, 

f^c = UjIi/u , (1) 

where cUpi = 47rasT^(l + Nf/6)/3 the color plasmon frequency||, |^ and z/ is a collision 
frequency. Color conductivity arises as a generalization of Ohm's law: — '^cEf and 
influences the evolution of a QGP through ohmic heating as pointed out in ref. [|TU|. Ohmic 



heating also couples mini-jets and beam jets in ultra-relativistic nuclear collisions as shown 
in ref. [0. The main problem in this connection is to estimate the collision frequency, v. 
In QED plasmas, v = 1/tp, where tp is the momentum relaxation time. In QCD, transport 
theory estimates for tp including color Debye and dynamical screening give|^ 

tp^iAa!\n{l/a,)T)-' . (2) 



This time scale was used to estimate the magnitude of effects of ohmic heating in ref. [ITJ] . We 
show below, however, that t^ instead of tp controls the color conductivity, and that tc ^ tp 
in the perturbative limit. In fact, we find that tc is closely related to the (divergent) gluon 
damping rate and that color diffusion strongly damps collective plasmon modes. We therefore 
find that non-abelian plasmas are surprisingly poor color conductors and that collective color 
modes are over-damped due to rapid color diffusion. 

To derive tc we follow here ref.|T2|, |13[ and consider the motion of a classical colored 
particle moving along a world line x^{t) with momentum p^{t) = mv^{T) and = dx'^/d.T, 
that couples with an A^^ — 1 component color charge vector A'^(r) to a background Yang-Mills 
field, We follow this approach here because of its intuitive simplicity and because it 

leads to the same result that we find starting from the kinetic equations for the quark and 
gluon Wigner density matrices [|l| including the Lenard-Balescu type collision terms derived 
in |1l4| , |15| . While the derivation based on the later approach is more rigorous, it is also less 



transparent physically and therefore deferred to a subsequent paper [Q. 

Consider the motion of a heavy quark in the d dimensional representation of SU{N) 
with second order Casimir, C2. The classical color charge vector is defined through the 
expectation value of the color current via 

rix) = g{i,{x)reij{x)) = g I dTv^{T)5\x - x{r))X^{r) , (3) 

where = if°''^H'^ and = C2I. Similarly, the classical momentum is related to the 

expectation value of the matter energy-momentum tensor 



T^'''{x) = i{^YD''{x)^{x)) = mj dTV^{T)v''{T)S\x - x(r)) , (4) 
where D'^{x) = d'^ + igt°'A'^^{x) The equations of motion follow from the conservation laws 

d.T^" = gf.Fr , d,j>^^ = gr'^Alr ■ (5) 
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Substituting Eqs. (^,^) into the left and right hand side, we recover the equations derived in 
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^ = gVuFrA-'' , (6) 
= gr'^^A'^A^ (7) 

where the fields are evaluated at x^(r). Since under a gauge transformation U{x) = 
exp{i6°'{x)t°'), ip Uip, the color charge must transform as t°A'^(r) U {x{T))t°' A"' {t)U~^ {x{t)) , 
to insure the gauge covariance of eqs.(H,|^). (Note that md/dr = p^d^ along the world line.) 

We study the time dependence of the color charge averaging over an ensemble of back- 
ground field configurations. We assume that the ensemble characterizes a color neutral 
medium in which the are random fields such that the ensemble averaged potentials van- 
ish in a suitable gauge. The ensemble average of Eq. (|^) gives 

^mr)) = gr''^{[vAtA%r)) , (8) 

where for shorthand we write [vA]'^ = {p^{r)/m)A^^{x{T)). Integrating (formally) Eq. (^ 
and substituting the result back into Eq. (|[), yields 

^(A'^(r)) = g'r'j^'^ J^^ dT'{[pAt[pA]t,A%T')) (9) 

Our first main physical assumption is that the time scale of field fiuctuations is short com- 
pared to the time scale of variations of the color orientation. Formally, this is suggested 
from eq.(^ if g <^ 1 because the rate of change of the color charge is proportional to g^. 
We therefore assume a stochastic (random phase) ansatz||l^ for the expectation value of the 
product 

(M]tM]f,A^(r'))^(M]tM^,)(A^(r)) . (10) 
With this ansatz, eq.(||) reduces to 

:(A"(r)) ^ {/r''T'^^|^^rfr'(M]tM]^,)}(A^(r)) = -d-(r)(A^(r)) . (11) 



d_ 

dr 



This defines the color diffusion tensor, (i'^*^(r), whose physical interpretation will become 
clear below. Note that we cannot use the v^A^^ = gauge in Eq. (|ll|) because ^^(t), via 
eq.(|^), is different for every member of the ensemble. 

Our second main approximation is to work in the eikonal limit. This is equivalent to the 
"hard thermal loop" approximation in We thus assume that our test parton has a high 
initial four momentum and that energy loss via eq.(^ is small. This also follows formally 
from eq.(^ in the perturbative <C 1 limit. In this case, we can factor the four velocity out 
of the ensemble average and approximate 

{[vAtivA]^,) ^ v'^v''{A'^{x{T))At{x{T'))) ^ v^v''Clt{x{T), x{t')) . (12) 
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Next we note that in this eikonal hmit x'^(r) — x^(r') = (r — t')v^^ and thus for a homoge- 
neous, color neutral ensemble, the correlation function C can expressed as 



Cj^,(a;(r),a;(r')) = 5m I {dk)e-'^^-^^^-^'^C^,{k,u) , (13) 



where {dk) = (i^fc/(27r)^, is the four velocity of the ensemble rest frame. The correlation 
function C measures the spontaneous fluctuations of the background field. In the classical 
limit those fluctuations can be calculated using kinetic theory 0, |T^. For a system in 



thermal equilibrium, the -C T spontaneous fluctuations can also be related to the response 



function (retarded commutator) via the fluctuation-dissipation theorem [|T9| as 



2T 

C"'ik,u)^-—lmD^^^ik,u,T) , (14) 

where 

D^j^ {k,u,T) = i j d^xe''="^(xo)Tr(p[A^(a;),A'^(0)]) 

is the thermal averaged retarded commutator that arises in linear response theory. The 
statistical operator is given by p = exp(— P^m^/T) for a thermal ensemble with four ve- 
locity u. Both methods give the same result in the high temperature limit in the 1-loop 
approximation]^. The 1-loop result obtained in ref.p, ^ is 

(15) 

where a is a gauge parameter. This decomposition utilizes the longitudinal and transverse 
projectors, Q^^ = u^u^/u^ and P^^ = g^^ - Q^^, as given in terms of g^^ = g^^ - k^k^/k"^ 
and = g^vU^ ■ Furthermore, the longitudinal and transverse polarization functions are 
related to the gluon self energy W{k) through Hl = H'^'^Q^^ and Ut = 
With (0,^), the color diffusion tensor can be expressed as 



d^^ir) = -2g'NT6ae J ^^ImivD^v) dr'e''^'^^^^-^'^ , (16) 

where the color factor for SU{N) follows from y^^cjefec _ ]\f§^^^ As r — oo, the proper time 
integral reduces to n6{kv) as can be seen using lmDji{—k) = —lmDji{k), to replace the 
limits of integration by — r to r. We can therefore define the color diffusion coefficient, dc, 
and the corresponding color relaxation time in the plasma rest frame, tc, via 

lim d^'ir) = 6aeNdc = 5ae{vu)/tc . (17) 

For ultra-relativistic partons, the color relaxation rate reduces to 

1 _ o_ 2..^ / W k^ Jkv\^^f 1 1 



-27rg^NT J 



(5 — Im 



tc "J {ku){kuy-k'^ \vu) VP -Hi P-Ht 

It is important to note that the dependence on the gauge parameter a dropped out because 
the delta function constrains kv = 0. However, in the semi-classical limit (a; -C T) the 
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gauge invariance of tc holds even if kv ^ 0, as can be shown via kinetic theory |T6[, since 
fluctuations are driven by gauge invariant current- current correlations (polarization tensor). 
The gauge invariance of this result is related to the gauge invariance of the damping rate in 



hard thermal loop calculations in ref.[^ and depends essentially on the eikonal limit. For 
momenta ~ 3T, the eikonal approximation applies perturbatively because the energy loss 
per interaction, ~ g'^T, is relatively small if g <^ 1. 

Given the solution of eq. ([TTp in the plasma rest frame is 

(A"(t)) = e-*/*=(A"(0)) . (19) 
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Note, that while the ensemble averaged color of the parton vanishes rapidly for t > t, 
the equations of motion conserve the magnitude of the color vector, i.e., (A°'(r)A"(r)) = 
constant. 

We can compare the color relaxation rate, 1/tc, to the quark or gluon damping rates 
notingpl that the longitudinal and transverse polarization functions are related to the color 
dielectric functions via fc^ — 11^ = k'^ei and k^ — Ut = u'^eT — k'^- The analytic expressions for 
eiik) and erik) are recorded for example in [|l], @]. In terms of these color dielectric functions 

ijj=k cos 



tc 27r2 J cjP \eL{k) u^rik) - k^ 



This is the same integral that occurs in the quark damping rate, 7g, as discussed in ref.0. 
In fact we find that it corresponds exactly to the leading log approximation to the gluon 
damping rate computed in This coincidence of results based on classical kinetic theory 
and 1-loop high T QCD is a general result [0] in all problems where the classical uj <^T and 
k <^T modes of the system are dominant. 

Unfortunately, as is well known [pO[|, the transverse contribution to that damping rate is 
logarithmically divergent. Dynamical screening is sufficient to regulate the Coulomb diver- 
gence in the expression for the momentum relaxation rate[^ but not the damping rate. The 
longitudinal part is finite and can thus be neglected. Several attempts to deal with 
this problem have been proposed. One is to introduce a nonperturbative magnetic mass, 
fTT'M ~ g'^T, as in ref. Another is to introduce damping self consistently by adding an 



imaginary part to the bare fermion propagator p2| . However, no satisfactory solution is yet 
available . 

In our kinetic formulation the source of the problem may be traced back to our first 
assumption that the field fiuctuation time is short compared to the color diffusion time. What 
we learned a postiori was that fluctuations of the transverse magnetic flelds are long ranged 
in time and that quasi-static unscreened magnetic flelds lead to the divergent color relaxation 
rate. This suggests that the change of the color moment with time should taken into account 
within the integrand in eq.(|TTp. Such a "memory" effect could damp the contribution from 
early times with {vu){t — t') > t^. Inserting a factor exp{—[vu){T — T')/tc) within the integral 
over t' would smear out the delta function, 6{kv), appearing in eq.(|TB|) into a Lorentzian 
form oc dc/^ikv)"^ + d^). The resulting self consistent equation for dc would lead to a flnite 
result analogous to the method considered by Altherr et al. p2 . 



We will however follow Pisarski[EO| and regulate the divergence by introducing a magnetic 
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mass cutoff. In the region u <^ k, we therefore approximate 

lm( ^ ]^ ^ (21) 

\u;hT{k)-ky {k^ + mlj^ + iau/ky ' ^ ^ 

where a = SirUpi/A ^ rn^. In the m^/ = limit this coincides with the one-loop form used 
for example in [^]. Note that the dynamic screening factor is properly included above. 
With eq.(^) we can integrate the transverse part to obtain 

1 ttgNT f dkk^dii auk'^ , , , , , 

" —J (P + ml,f + (a,f '' "-"^ ^osim.M , (22) 

where the final form is numerically accurate already for ttlm < uie- The finite longitudinal 
part adds about 0.5 to the log which is well within the uncertainties associated with inserting 
a magnetic mass in (^I]) by hand. 

To show that dc actually corresponds to a diffusion coefficient in kinetic theory, we con- 
sider next the classical phase space density A) of an ensemble of colored particles 
defined by 

Q(x,p,A) = j dT'Y,5{x-x,{T'))5{p-p,{T'))5{K-K{r')) , (23) 

i 

where pf and A" obey Wong equations (|^) coupled via a self consistent field. The phase 
space density obeys the (Heinz) transport equation [1^ : 

M + OT'^i^;.A"5p^ + igv'AlL^)Q{x, p, A) = , (24) 

where 

are generators of SU{N)/Z^ obeying the commutation relations [L"", L^] = if"'^'^L'^. eq Note 
that ( p^ just non-abelian Louiville equation expressing that the parton density is constant 
along a characteristic, (x(r),p(r), A(r)), satisfying the Wong equation in the self consistent 
field, i.e.. 



d 

dr 



Q(x(r),p(r), A(r)) = - | dr'^ ^(^M " ^^(^')) " " " = • (26) 



The ensemble average of eq. ( PB]) can be expressed as 
d^ 



m—{Q{x{T),p{T),A{r))) = -ig{jf5AimQ) - g{p'5F;^K'^d>;5Q) , (27) 



where we decompose Q = (Q) + 6Q and A = {A) + 6A, and the characteristic trajectory 
satisfies 

dr^^ dr)^ d K"" 

'^^=P' ' ^ = ^7^.(^r)A" , -^=gr'%^{Al)A^. (28) 



To proceed, we assume average color neutrality, {A) = {F) = 0, and neglect the momentum 
degradation part of the kinetic equation to focus exclusively on color diffusion. Therefore, 
the mean phase space density evolves according to 

m-^{Q{x{T),p,A)) = -zg{v^6A';^L'^6Q) . (29) 

The kinetic equation for the fluctuating part of the phase space density is obtained by 
subtracting eq.(p7D from .(pG]): 



d 

dr 



6Q{x{t),p{t),A{t)) = -^g[v6Ar^L'^{Q) - gv''6F;;^A^d^{Q) , (30) 



where in this case 

= gv^HFn + 5Fn^- , = gr^'-HAl) + 54)A'= . (31) 



Taking into account our assumption (A) = (F) = 0, neglecting the induced fluctuations 
caused by the gradient in momentum, and retaining only terms linear in fluctuating quantities 
Eq. (|30|) reduces to 

^SQ{x{t),p,A) = -zg[vSAr,L'^{Q) , (32) 

This linearized approximation implicitely assumes the eikonal limit since then p is constant 
and also A is constant to that order. Thus both the assumptions made in deriving the color 
diffusion coefficient previously correspond to this linearization of the kinetic equations for 
fluctuations. Integrating eq.(|32|) and substituting the induced fluctuation into eq.(p^) leads 
to 



d{Q) 2 



g' dT'{[v6AUv6At,)L^L''{Q{x{T'),p,A)) . (33) 



dr 

Assuming that fluctuation correlation time is short compared to the color diffusion time, 
{Q{x{t'),p, A)) ^ {Q{x{t),p, A)) allowing us to factor Q out of the integrand. In this way, 
we obtain for homogeneous and color neutral ensemble the following color diffusion equation 

^(g(x(r),p. A)) = -d^L'^L^Q) , (34) 



where color diffusion coefficient dc is given by 

4 = (^V^-r rfr'C^,(x(0),x(r')) , (35) 



where C is given by eqs.(|T2|-|I3|), is thus the same color diffusion coefficient that we derived 
via eqs.(|16|- 0)). Note that Eq. (|3^) can be written in the form of Fokker-Planck equation 
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where color diffusion coefficient d""^ is given by 



Equation (p4D therefore describes diffusion in color space. Note that operator C2 = L"-L"- 
just is quadratic Casimir operator of SU{N)/Z'^ . For a more detailed derivation of collision 
terms including momentum space diffusion we refer to fl^, |15|, p!6| . 

To illustrate color diffusion, it is instructive to consider SU{2), where = —le^^^K^d/dA.^ 
is an angular momentum operator in color (iso-spin) space and C2 = l{l + !)• Since the 
magnitude of A is fixed, any initial distribution in color space can be expanded in the 
spherical harmonics as 

(Q(A,r = O) = 5:Q,„Fr(^,0) , (38) 

where 9, are polar coordinates in iso-spin space. In this case, the proper time evolution is 
simply 

(Q(A, r)) = ci,mYr{0, 0)e-^-'('+i) . (39) 

1,'m 

From Eq. ([39| ) it follows that any non-isotropic distribution in color space evolves to a uniform 
distribution, 

J[im(Q(A,r)) = co/v/47r, (40) 

on a proper time scale ~ l/2dc (for = 2). The same picture holds qualitatively for 
the general SU{N) case, where any distribution evolves to uniform distribution which is 
determined by the zero mode of quadratic Casimir operator. This color equilibrium state is 
achieved after a proper Tc = 1/Ndc- Since the time in the plasma rest frame is t = {vu)t 
and Ndc = {vu)/tc, color equilibrium in the plasma rest frame is achieved on the time scale 

Finally, to calculate the color conductivity coefficient, we must test the linear response 
of the system to a weak external field. In an external filed, F^x, the kinetic equation for the 
mean density must be supplemented by a Vlasov term, gpFex-^dp{Q) , on the left hand side 
of Eqs.(|^,0). To calculate the induced color current, we multiply that modified Eg . (p4|) 
with A" and integrate over the color charge with the measure pTSf 



dAoc II rfA,5(A"A, - q2)6{A''A'A'dabc - gs) ■ (41) 
1=1 

The proportionality constant is chosen such that [T^ / dA = 1, / dAA"- = 0, and / (iAA'^A* = 
C2(5"^ The char ges are defined through the trace of the Casimir operators: q2 = Tr(t'^t"), 
qs = dabJTTitHH^). For quarks ga = (A^^ - l)/2, gs = {N^ - 4)(Ar2 - I) /AN and ca = 1/2. 
For gluons = N{N'^ - 1), q^ = 0, Ca = A^. 

This leads to the following linearized kinetic equation for the adjoint color moment den- 
sity, g"(a;,p), 

v^d^q'^ + C2gv^[F,x]l,d;{Q), = -Nd,q- . (42) 
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where 



g«(p,x)= dAA^{Q{x,p,A)) , (43) 



For a homogeneous, color neutral ensemble, {Q{p,x,A)) = {Q{p))o, and thus since / dAA"- = 
the color moment density vanishes. However, in a weak external field there is a distortion 
of the color density linear in F^x- The factor on the right hand side of ( ^2]) follows 
from / dAA'^LF'Q = Nq"", which can be verified using [L",A''] = if°'^'^A^ and the property 
that / dAL°'F{A) = 0, which follows upon integration by parts and noting the symmetry 
properties of f"'^'^ and d"''"^. We note also that in contrast to (42) the kinetic equation for the 



color singlet density, q{x,p) = J dA{Q), is independent of the color diffusion coefficient. 
Solving (1^) for the Fourier transform, q°'{k,p), the induced color current is given by 

r{k) =g I dpp^^q^ip, k) = a^"^(A;)[Fe,(A;)]^^ , (44) 



where the conductivity tensor for a homogeneous color neutral plasma of quarks, antiquarks, 
and gluons is 0, 

a^^(t) ^ «7^P _.(^,f;;^);,/ ?{^(/.M + fM) + mm . (45) 

where, fi{p) are the momentum space distributions including spin degeneracy for the partons. 
For example, in equilibrium (p) = (2/(27r)^)20(po)^(p^)/(e^"''^— !)• For an isotropic plasma 
this tensor reduces to a^"'^ = a^'^{k)uf^ . Finally, from Eq. (^5]) it follows that 0] the static 
[k = 0) color conductivity is cr^'^(O) = cTcg'^'^, where 

cTc = t,uli . (46) 

The important difference between our derivation and that in ref. @ is that the collision rate 
is derived here directly from the kinetic equation for color fiuctuations rather than parame- 
terized via the relaxation time approximation. We see in particular that the color diffusion 
time rather than momentum degradation time controls the magnitude of the conductivity. 
In addition, the same color diffusion rate applies for both quarks and gluons unlike in the 
general ansatz of Inserting the expression for tc, we thus find the surprising result stated 
in the introduction, i.e., that 

a, ^2T/log(ms/mM) • (47) 

The appearance of the non-perturbative and non-classical magnetic mass in this final ex- 
pression suggests that the classical {k —>■ 0) or high T limit of a QGP is in fact be very 
different from QED plasmas because of color diffusion. Another indication of this is that 
the long wavelength [k <C T) color plasmon mode, obtained from the dispersion relation 
/c^ = Hiik) = —iQ ^v{'^^°"' — (J'^'^°')ka, are strongly damped with a rate proportional to 1/tc 
as can be seen from the general discussion in ref. . Thus, the anomalously large damping 
rate of hard thermal partons propagates down to the soft collective modes. Probably not only 
the plasmon but also the plasmino are over-damped due to color diffusion. In this respect, 
further study of color diffusion and screening of transverse interactions are very important 
in the development of QGP transport theory. In practical applications, especially to nuclear 
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collisions, of course all time integrals are restricted to a finite time, which automatically reg- 
ulates the formal static divergences. The problem in that connection is to determine which 
effect, finite time or nonperturbative magnetic screening, restricts the magnitude of color 
diffusion the most. 

Acknowledgments: We thank R. Pisarski, S. Gavin, and T.Ludlam for fruitful discussions. 
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